The breakdown of the zeroth law of thermodynamics and the definition of 
temperature in small quantum systems 
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We study two small quantum systems coupled to the same reservoir which is in thermal equi- 
librium. By studying the particle density and the energy density in the two systems before and 
after they contact each other, we find that the two systems are not in thermal equilibrium with 
each other. Our result shows that the zeroth law of thermodynamics is broken in small quantum 
systems at low temperatures. Therefore, the traditional way of defining temperature fails due to 
the breakdown of the transitive relation of thermal equilibrium. Then we show a different way of 
defining temperature by attaching an auxiliary site, which plays the role of a thermometer, to the 
small quantum system. 



I. INTRODUCTION 

Thermodynamics, developed in the 19th century, con- 
tains four fundamental laws. Among them the zeroth 
law of thermodynamics states the transitive relation of 
thermal equilibrium and is the prerequisite of the exis- 
tence of temperature. For more than one hundred years, 
it is believed to be a fundamental law of nature and valid 
throughout various systems, either classical or quantum. 

The zeroth law of thermodynamics states that two sys- 
tems in thermal equilibrium with a third one are also in 
equilibrium with each other. In other words, the ther- 
mal equilibrium is a transitive relation, so that one can 
define a physical quantity called temperature for a set of 
systems which are in equilibrium with each other. This 
statement can be checked theoretically by studying the 
dynamics of three systems coupled together. In detail, 
we suppose two systems coupled to a large thermal reser- 
voir at the temperature T. At the initial time, the two 
systems are isolated to each other and are both in equi- 
librium with the reservoir. According to the zeroth law of 
thermodynamics, they must also be in equilibrium with 
each other. This can be checked by switching on a cou- 
pling between the two systems. If they have been in 
equilibrium, there should be no exchange of energy or 
particles between them after the coupling is switched on. 
Conversely, if switching on the coupling causes the ex- 
change of energy or particles, which can be shown by the 
changing of energy density or particle density in the two 
systems, the zeroth law breaks down. 

In this paper, we check the zeroth law in small quan- 
tum systems which obey the quantum laws. In recent 
years, with the development of the nanotechnology in- 
tense efforts have been devoted to study the thermaliza- 
tion process in different small quantum systems. Studies 
show that an open small quantum system may not evolve 
into an equilibrium state with a canonical density matrix 
as the macroscopic systems [H4ri| . This result challenges 
the method of measuring the temperature. Because the 



standard way of measuring the temperature of a system 
is to attach a thermometer to it. The thermometer will 
then evolve into an equilibrium state with the tempera- 
ture as same as the system, which can then be readout. 
However, if the thermometer does not thermalize, the 
readout of it will depend upon how it is coupled to the 
system and then lose the meaning of temperature. 

In this paper we will show that not only the method 
of measuring temperature but also the definition of tem- 
perature is challenged in small quantum systems. The 
traditional definition of temperature is based on the ze- 
roth law of thermodynamics, which will be shown not to 
stand up in small quantum systems. We will introduce a 
different way of defining the temperature, which circum- 
vent the zeroth law by attaching an auxiliary system to 
the small quantum system. 

The paper is organized as follows. In Sec. II, we intro- 
duce the main model in this paper, an electron reservoir 
coupled to two small systems, and the method of solving 
it. In Sec. Ill, we discuss the breakdown of the zeroth law. 
A method of defining the temperature in small quantum 
systems is introduced in Sec. IV. In Sev. V, we extend 
our discussion in a model where the number of particles 
is conserved. Sec. VI is a short summary. 



II. MODEL AND METHOD 

We employ two chains consisting of n sites respectively 
to model the two small quantum systems, which are cou- 
pled to an electron reservoir at the temperature T. At 
initial time the two chains are isolated to each other. To 
see if they have been in thermal equilibrium with each 
other, we switch on a coupling between the end sites of 
the two chains. The Hamiltonian describing the chains 
and the reservoir is expressed as 

H = ^ tkc\ a c ka + yfoa4q g <*lg + h.c.) 



n-1 



*Electronic address: wangpei@zjut.edu.cn 



+ e a L<T 5 a J <r+ 9a(4iaCa,3+l,a + h.C.) 

(y,(J.j — l Oi(T,j — l 

+d(c[ n<T c Rnry + h.c), (1) 



2 



where a = L,R denotes the two chains, a =f, I the spin, 
Cfc CT the annihilation operator of the electron in the reser- 
voir, c a j a the annihilation operator at site j in the chain 
a, h a the coupling between the chain a and the reservoir, 
g a the coupling between two neighbor sites in the chain 
a and d the coupling between two chains. 

The initial correlation is absent in this model. So 
the Keldysh Green's function approach [l8|, [l9[ can be 
employed for calculating the electron density in the two 
chains before and after the coupling d is switched on. 

The spin-dependent electron density at arbitrary site 
to in the chain a can be expressed by the lesser Green's 
function as 
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The Gumam is the element of a 2n x 2n Green's function 
matrix, which according to the Keldysh formalism can 
be related to the retarded and advanced Green's function 
matrix by 



G < 



G r T, < G a . 



(3) 



In this formula, the lesser self-energy matrix contains 
only four non-zero elements, which are listed below: 



{ 



e< >L1 = 2ir L /H 

S L1„R1 = S fl.l,Ll = 2TW(w) 

1 



(4) 



Here f(uj) — g denotes the Fermi distribution 

e k B T + 2 

function, r a = /07r/i„ and = p-KhjJi R denote the 

level broadening at the edge sites due to the coupling to 
the reservoir and p is the density of states in the reservoir. 
Substituting Eq. @] into Eq. [3J we obtain 



Gam,am ~ 2*/ ( I G a m,Ll I + \Gam,Rl\ ^1 

+2T LR^G.(G r am L1 G r ^ m R1 )) , 



(5) 



where we use the relation G° = G r ^ . Then the electron 
density can be written as 
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-2T LR Ke{G r am Ll G r ^ m Rl )) 



(6) 



All the 2n x 2n retarded Green's functions are the el- 
ements of a matrix G r , which satisfies the Dyson's equa- 
tion: 



G r (uj) = G 0r (uj) + G 0r {uj)T, r {uj)G r (Lj). 



(J) 



Here G 0v (lu) is the free Green's function matrix as h a 
g a = d = 0, and its elements are 



G a r j , a > f (u) = 6 jtj ,5 a „ 
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uj — e + ir\ 



(8) 



where 77 is an infinitesimal number. And S r (w) is the 
self-energy matrix, whose elements come from the kinetic 
energy of electrons hopping between neighbor sites and 
the level broadening at the edge sites. The subdiagonal 
and superdiagonal elements of Y7 (uS) are 

^aj,aj + l = ^aj+l.aj = 9a- (9) 

The coupling between the left and right chain gives 

^Ln,Rn = ^Rn.Ln ~ d- (10) 

And the level broadening at the edge sites gives 

r E^ 1>al (a;) = — iT a , . 

Solving the Dyson's equation, we obtain the retarded 
Green's functions: 



G r = (G 



Or-l 



57)- 



(12) 



Finally, the electron density is got by substituting Eq. fTS] 
into Eq. [5] 

In the above formalism, the electron-electron interac- 
tion is not considered. The interaction can be included by 
adding a Hubbard term H v = U J2 a j ^ajt^t^aji^aji 
to the Hamiltonian. In the self-consistent mean field ap- 
proximation, the interacting model can be easily solved 
by replacing the on-site potential e by the e QJ - = e + U n a j 
in the above formalism and then calculating the electron 
density n a j self-consistently. This has been proved to be 
a good approximation as the interaction U is small. 



III. THE BREAKDOWN OF THE ZEROTH 
LAW OF THERMODYNAMICS IN SMALL 
QUANTUM SYSTEMS 

First we consider the special case that each chain con- 
tains only a single site and the interaction [7 = 0. Now 
the chain L and R are also called the site L and R re- 
spectively. The spin-dependent electron density at the 
site L is 



du- 



dVr^Y 



and that at the site R is 



n R = — l du- 
ir 



^(w - 



dJT, 



\M 2 



/M, (13) 



-/M. (14) 



where A = (cj - e + iT L )(uj - e + iT R ) - {iT LR - d) 2 . 

The electron density at the sites depends strongly upon 
the coupling d as the level broadening of the two sites are 
very different. We set the Fermi energy of the reservoir 
to be the energy zero. In the case of the on-site energy 
level e < and T L >• |e| ^> T R , the electron density at 
the strongly-coupled site (the site L) is = 0.5 before 
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FIG. 1: The electron density at the weakly-coupled site R as 
a function of the coupling d at different temperatures. The 
Fermi energy is set to be the energy zero. The level broad- 
ening at the two sites is Tl = 10r.L.R and T_r = O.ITlr re- 
spectively, and the energy level is e = — Tlr- The black circle 
denotes the electron density when the two sites are isolated 
to each other (d = 0). 
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FIG. 2: The electron density at the weakly-coupled site R 
as a function of the coupling d at different temperatures and 
the energy levels e. The level broadening is Tl = WTlr and 
r_R = O.ITlr- The energy level e takes two possible values: 
and Tlr- And the temperature is chosen as T — OK, 0.5Tlr 
and Tlr- The black circle denotes the electron density when 
the two sites are isolated to each other (d = 0). 



the two sites contact each other. And the electron den- 
sity at the weakly-coupled site (the site R) is ur = 0, 
since it is screened. Switching on a coupling between 
the two sites has little effect on the electron density at 
the strongly-coupled site. But the electron density at 
the weakly-coupled one will change abruptly into a finite 
value depending upon the temperature. As d increasing 
the electron density at both sites will go towards 0.5 (see 
Fig. 1). 

We see that the electrons at the two sites will redis- 
tribute themselves after the sites contact each other. So 
the two sites which are in equilibrium with a reservoir are 
not in thermal equilibrium with each other. The transi- 
tive relation of thermal equilibrium is broken in this case. 
This picture will not change as e > (see Fig. 2). In a 
conclusion, the electron density at the weakly-coupled 
site will change significantly after the two sites contact 
each other whatever the temperature and the on-site en- 
ergy level are. The two sites are not in equilibrium as 
their level boadening is much different from each other. 

The breakdown of the zeroth law is because the quan- 
tum correlation between the two sites is sensitive to the 
coupling strength d. But as the chain length increasing, 
the coupling at the boundary will become unimportant. 
And then the zeroth law should recover. This is verified 
by the numerical calculations. 

We study two chains which are coupled to the reser- 
voir in the very different strength, i.e., Tl ~> Tr- The 
averaged electron density in the weakly-coupled chain 
as a function of d is shown in Fig. 3 for different chain 
length. We find that the electron density changes with 
the coupling d. But the change becomes small as the 
chain length n increasing. 

And at the finite temperature, the thermal fluctuation 
will suppress the influence of the coupling between the 
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FIG. 3: The averaged on-site electron density in the weakly- 
coupled chain as a function of d at different temperatures. 
The temperature is chosen asT = OK, O.ITlr and Tlr- The 
length of the chains varies from 6 to 20. The chain L is 
strongly-coupled to the reservoir with Tl — Ql — \0Tlr, 
while the chain R is weakly-coupled to the reservoir with Tr = 
gR — 0.ir_Li?. The energy level is set to be e = — Tlr- As 
the chain length increasing, the dependence of the electron 
density on d is reduced. 



two chains. This is clear in the case of both n = 1 and 
n > 1 (see Fig. 1, 2 and 3). The breakdown of the zeroth 
law is only significant at low temperatures. 

We study the effect of the interaction U. As shown 
in Fig. 4, the interaction changes the electron density. 
However, what we are interested in is whether switching 
on a coupling d will change the electron density. This is 
found to be independent to the interaction. And we also 
find that the hopping g L and gn in the two chains is not 
important in our picture. 
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FIG. 4: The averaged on-site electron density in the weakly- 
coupled chain as a function of d at different U, gL and ga- 
The temperature is chosen as T = Tlr- The length of the 
chains is n = 6. And the coupling between the chains and the 
reservoir is Yl = WTlr and Tr — 0.ir_L_R respectively. The 
energy level is set to be e = — Tlr. 



In a word, the transitive relation of thermal equilib- 
rium breaks down in short chains at low temperatures. 
This indicates that the traditional method of defining 
temperature fails for small quantum systems at low tem- 
peratures. It is necessary to find a new definition of tem- 
perature in this situation as the correlation at the bound- 
ary cannot be neglected. As will be shown next, such a 
definition can be obtained by attaching an auxiliary site 
to the system that we want to measure. We also find 
that the zeroth law is recovered when we increase either 
the temperature or the system size. So the zeroth law is 
always a good approximation in everyday's life. 



IV. THE LOCAL TEMPERATURE IN A SMALL 
QUANTUM SYSTEM 

Due to the breakdown of the zeroth law in small quan- 
tum systems at low temperatures, we introduce an alter- 
nate way of defining the temperature, which is inspired 
by the process of using a mercury thermometer to mea- 
sure the temperature of an object. The measuring pro- 
cess is to attach the thermometer to the target and wait 
for a long time, so that an equilibrium between the two 
systems is established. The temperature of the target is 
as same as that of the thermometer, which can then be 
obtained from the relation between the temperature and 
the volume of the mercury. 

The condition that this procedure works is that the 
thermometer after attached to the target should be in a 
canonical state. This can be realized when the coupling 
between the thermometer and the target is infinitesimal. 
An example showing this is in the Anderson impurity 
model without interaction, which describes an impurity 
site coupled to a Fermi sea at the temperature T and the 
chemical potential fi. The impurity site plays the role of 



the thermometer. As is well known, the electron density 
at the impurity site goes to the Fermi function 



1 



'(«d-M) 



1 



(15) 



in the limit as the coupling between the site and the 
Fermi sea goes to zero. Here denotes the energy level 
of the impurity site. This fact indicates a method of 
defining the temperature of a small quantum system at 
low temperatures. 

We define the local temperature of a quantum system 
in thermal equilibrium by attaching an auxiliary site to 
the system. The auxiliary site plays the role of the ther- 
mometer. The coupling between the auxiliary site and 
the quantum system is switched on at the time t = 0. 
As the coupling goes to zero, the electron density at the 
auxiliary site as a function of its energy level must be the 
Fermi function when the time goes to infinity. According 
to Eq. [TS1 the electron density n^ej) satisfies 



In [ — 



1 



1 



knT 



(16) 



Then the temperature T and the chemical potential fi 
of the quantum system can be determined by plotting 
\n{l /rid — 1) with respect to e^. The slope of this linear 
function is just the inverse of the temperature * 

The temperature that we obtain in this way is me local 
temperature, which depends upon the position where the 
auxiliary site is coupled. For a system in thermal equi- 
librium, the local temperature is the same everywhere, 
equaling to the global temperature. 

Now let us study the local temperature of two sites cou- 
pled to an electron reservoir in thermal equilibrium with 
the different level broadening, in which we have shown 
that the zeroth law of thermodynamics is broken. 

The Hamiltonian describing the two sites, the reservoir 
and the auxiliary site is similar to Eq. [1] i.e., 



ka Oi'ka 



a' a <7 

+g'^2(dlc aa + h.c), (17) 

cr 

where d a denotes the annihilation operator at the auxil- 
iary site, ed the energy level of the auxiliary site and g' the 
coupling between the auxiliary site and the site a whose 
temperature is planned to be measured. For the given a 
(a = L, R), the electron density at the auxiliary site as a 
function of e d can be calculated by the Keldysh Green's 
function method. We find that the spin-dependent elec- 
tron density can be expressed as 



1 

n d = - 

IT 



duf(\G 



d.L 



2 r L 



\G dtR \ T R 



+2T LR Re(Gl L G d : R )) , 



(18) 



5 



J3 




FIG. 5: The spin-dependent electron density rid at the aux- 
iliary site as a function of its energy level. The black circle 
and the white square denote the electron density as the aux- 
iliary site is coupled to the site L and R respectively. The 
site L and R is coupled to a Fermi sea in equilibrium with 
the level broadening of Tl ~ WTlr and Tr = O.ITlr re- 
spectively, and the energy level of the two sites is e = — Tlr- 
The solid line denotes the Fermi function at the correspond- 
ing temperatures. We choose three different temperatures. In 
all the situations the function n^ed) fits well with the Fermi 
function. 



where G r d L and G r d R denote the retarded Green's func- 
tions between the auxiliary site and the site L and R re- 
spectively. These retarded Green's functions can then be 
obtained by solving the Dyson's equation. The electron 
density as a function of the energy level at the auxiliary 
site coupled to either the site L or the site R is shown in 
Fig. 5. We find that as g' — > the data does fit well with 
the Fermi function at the temperature of the reservoir. 
This means that our method of defining the local tem- 
perature is appropriate in small quantum systems where 
the correlation between the system and its environment 
cannot be neglected. 

It is well known that the partial density matrix of a 
small quantum system in thermal equilibrium is not a 
canonical one due to the coupling between the system and 
the environment. Our result shows that it is possible to 
assign a temperature to a system even its density matrix 
is not a canonical one. 



V. THE DEFINITION OF TEMPERATURE IN 
A SYSTEM WHERE THE NUMBER OF 
PARTICLES IS CONSERVED 

In above, we discuss the breakdown of the zeroth law 
in the systems consisting of identical fermions, in which 
we define the local temperature and the local chemical 
potential simultaneously because the systems are open. 

In this section, we turn to a closed system, in which 
the number of particles is conserved. Then we could tem- 
perarily forget the chemical potential and focus on the 



temperature. In details, the systems contain a series of 
spins with nearest-neighbor couplings. The reservoir is 
simulated by an infinite spin chain. A local magnetic 
field is supposed to split the energy levels of the spin up 
and spin down states in the systems. The Hamiltonian 
of the systems and the reservoir is expressed as 

/-(n+2) oo \ 

H=9 £ + £ {sfs! +1 + sfs! +1 ) 

\ i= — oo i=(n+l)/ 

a—x : y a—x,y 
-2 n-1 

+ gL £ £ SiSi+l + SR £ £ Si&i+l 

i—1 a—x,y 



%— — n a—x,y 



-1 n 



+ d £ S^Sf + el £ +E )Sf, 



a—x,y 



\i— — n i—1 / 



(19) 



where the sites i with |i| > n are in the thermal reservoir. 
The site at the infinity is supposed to be coupled to the 
site at the minus infinity, so that the reservoir is as a 
whole. The sites i with — n < i < —1 are in the system 
L, and that with 1 < i < n are in the system R. The hh 
and ha denote the coupling between the reservoir and the 
system L and R respectively, and d denotes the coupling 
between the two systems. And e denotes the Zeeman 
splitting energy in the systems. 

The Hamiltonian can be re-expressed in terms of spin- 
less fermions through the Jordan- Wigner transfor- 
mation: 

/-(n+2) oo \ 

V i= — oo i=(n+l)/ 

+ ~Y ( d -(n+lf-n + h - c ) + (4cn+l + h.C.) 

-2 n-1 

+ £ (4^+1 + h - c ) + y £ ( 6 ^ i + i + h - : 

i— — n i—1 

+ ~($-A + h.c) +e( +£) e !c,- 



\i——n i—1/ 



(20) 



We immediately notice that this model is exactly equiva- 
lent to the model described in Eq.[T] except that in Eq.[20l 
the Hamiltonian of the reservoir is expressed in the real 
space while it is expressed in the momentum space in 
Eq. [1] And the spin degree of freedom is neglected in 
Eq. [20l since it is irrelevant as the electron-electron in- 
teraction is absent. The two end sites of the reservoir are 
at the positions of ± (n + 1). And the level broadening 

h 2 

at the edge sites of the system L and R are Tl = -g^ 

h 2 

and Tr = ^ respectively, due to the couplings to the 
reservoir. The local observable now becomes the spin 
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density, which is directly related to the particle density 
by (Sf) = - |. 

Since this model is equivalent to the one that we dis- 
cussed above, we get that the zeroth law is broken in this 
model at low temperatures. And the traditional way of 
defining the temperature fails for these closed systems. 
According to the above idea, the local temperature is de- 
fined by attaching an auxiliary spin at the position where 
we plan to measure the temperature. The auxiliary spin 
is a two-level system. A local magnetic field is applied to 
it so that its energy levels are splitted. Then the temper- 
ature of the auxiliary spin can be readout from the spin 
density as a function of the Zeeman splitting energy A. 
The Hamiltonian describing the auxiliary spin and the 
coupling between it and the system is expressed as 

H a = AS%+g' J2 SoS?- (21) 

a— x,y 

Without losing generality, we suppose that the auxiliary 
spin locates at the site and is coupled to the site 1 in the 
system R. After the Jordan- Wigner transformation, the 
model changes into the one described in Eq. [T7] except 
an unimportant interaction term. As g' — » 0, the spin 



density at the auxiliary site will be 

Obviously, the local temperature can be readout by plot- 
ting So as a function of A. 



VI. CONCLUSIONS 

In this paper, we show the breakdown of the zeroth 
law of thermodynamics in small quantum systems at low 
temperatures. The traditional method of defining the 
temperature, in which the temperature is a quantity as- 
signed to a set of systems in equilibrium with each other, 
fails in small quantum systems. Then we show a differ- 
ent way of defining the temperature. In this way we can 
obtain the local temperature of a system in equilibrium 
by attaching an auxiliary site to it. The temperature 
defined in this way is not an intrinsic quantity of the 
system's density matrix any more. 
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